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We study in detail a system of two interferometers aimed to the detection of extremely faint 
phase-fluctuations. This system can represent a breakthrough for detecting a faint correlated signal 
that would remain otherwise undetectable even using the most sensitive individual interferometric 
devices, that are limited by the shot noise. If the two interferometers experience identical phase- 
fluctuations, like the ones introduced by the so called “holographic noise”, this signal should emerge 
if their output signals are correlated, while the fluctuations due to shot noise and other independent 
contributions will vanish. We show how the injecting quantum light in the free ports of the interfer¬ 
ometers can reduce the photon noise of the system beyond the shot-noise, enhancing the resolution 
in the phase-correlation estimation. We analyze both the use of two-mode squeezed vacuum or 
twin-beam state (TWB) and of two independent squeezing states. Our results basically conhrms 
the benefit of using squeezed beams together with strong coherent beams in interferometry, even 
in this correlated case. However, mainly we concentrate on the possible use of TWB, discovering 
interesting and probably unexplored areas of application of bipartite entanglement and in particular 
the possibility of reaching in principle surprising uncertainty reduction. 

PACS numbers: 42.50.St, 42.25.Hz, 03.65.Ud, 04.60.-m 


I. INTRODUCTION 

The possibility of increasing the performances of interferometers by using quantum light represents one of the most 
interesting use of quantum states for overcoming classical limits of measurements. The first approach proposed to 
this aim is based on exploiting squeezed light for reducing the noise level in interferometers [T] and found recently 
application in gravitational waves detectors Hi. A second approach is considering the use of entanglement in phase 
estimation and, in particular, the possibility offered by the use of NOON states |1]. However, even if this approach 
presents a signihcant conceptual interest and could hnd very interesting applications in the future, nowadays the 
difficulty in producing high N entangled states and the fragility to noise and losses of these schemes strongly limits 
their real possible use. More recently, correlation in photon number in two-mode squeezed vacuum or twin-beam 
state (TWB) [5] has been demonstrated to be an important tool for beating shot noise H realising a first 

quantum protocol effectively robust against noise and losses [Tj- These results prompted to study the possibility of 
improving the so called “holometer” by exploiting quantum light, and in particular squeezed one or TWBs [5]. The 
bolometer is a double Michelson Interferometer (MI) addressed to detect the so called “holographic noise” (HN), 
namely a basic form of noise conjectured in quantum gravity theories that would derive from a non-commutativity 
of the spatial degrees of freedom at the Planck scale [9]- This noise, albeit very small, should be correlated when the 
two Mis are parallel, such to be in the respective light cones, while should be uncorrelated when one arm is rotated 
to be oriented in the opposite direction for the two Mis. The evident huge impact that the discovery of holographic 
noise, the first eventual evidence of quantum gravity effects, would present, [lOj motivates an accurate analysis of the 
possibility of improving the holometer performances. In this paper we detail and complete the analysis of Ref. [8] 
identifying operative situations where the use of quantum light would allow to greatly increase the performances of a 
double interferometer like the holometer. 

More in detail, in [8] we investigated an unusual but potentially powerful system consisting of two interferometers 
whose correlation of output ports signals is measured (see Fig0. This kind of double interferometric system can 
represent a breakthrough for detecting a faint correlated signal that would remain otherwise hidden even using the 
most sensitive individual interferometric devices, limited by the shot noise. On the other side if the two interferometers 
are in the experience identical fluctuation, this signal should emerge in a correlation measurement of their output, 
while the fluctuations due to shot noise and other independent contributions will vanish. 

The first experimental realization of this scheme using coherent beams (stabilized lasers), exactly for HN detec¬ 
tion, is already being implemented at Fermilab. Other applications can be envisaged such as for new generation of 
gravitational wave detectors. 

In [S] we have introduced a rigorous quantum model for describing the system. As opposed to standard phase 
measurement in a single interferometer, which involved first order expectation value of the output, in the double 
interferometric scheme the quantity under estimation is the covariance of the two outputs, which is a second order 
quantity, thus the associated uncertainty is a fourth order function. Notwithstanding this difference, we demonstrated 
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how the injection of quantum light at the input ports, which would remain unused in classical Holometer configuration, 
can boost the sensitivity of the device. In addition to the classical intense coherent beam, we considered both the use of 
independent squeezed beam (SQB) and correlated state such as the TWB. The ideal experiment described theoretically 
there, is however not suitable for a practical experimental implementation. In particular, the setting of the central 
phases the two interferometers exactly at = (f> 2 ,o = 0 which provides the optimal quantum enhancement, is indeed 
critical, because minimal deviations from this working regime completely compromise the advantages of the quantum 
strategy. Furthermore, the balanced readout configuration explored in SQB case would require simultaneously high 
dynamic range, fast and high resolution detectors that are not yet available. 

Here we present a framework in which a more complete and general study of the double interferometric system is 
provided, leading also to the depiction of a more experiment-oriented configuration of the system in terms of readout 
strategy and parameters choice. 

In Sec. [IT] we present the description of the scheme. In Sec. m we analyse the correlation properties at the output 
ports of the interferometer in the case of TWB, demonstrating that for a proper choice of the relative phase "0 between 
the coherent beam and the TWB, the energy transfer from the coherent beam due to the mixing allows the generation 
of either bright quantum correlation (^ = 7r/2) or anti-correlation (■i/i = 0) in the photon number, as witnessed from 
the Noise Reduction Factor (NRF) value below unity. Furthermore, we demonstrate that there are two regimes with 
different behaviour of the system: (A), in which a strict choice of the central phases of the interferometers makes 
the contribution to the coherent light at the dark port negligible, i.e. only quantum light is detected. In principle 
it leads to the optimal correlation (only bounded by losses), but turns out quite challenging in practice, since would 
require extremely precise stabilization of the interferometers, especially if the coherent beam is intense.(B), in which 
the major component of the signal at the output ports is due to the coherent beam contribution. This is a more 
common and realistic working condition. There is a sudden transition between these two regimes. 

In Sec. |IV| we describe in detail a model establishing the connection between a generic measurement operator 
(observable) and the estimation of the phase-covariance introduced by a correlated faint phase signal such as HN. 
Then, we focus on two specific quantum strategies: either on the use of TWB state and the measure of the photons 
number difference, of the use of two independent squeezed states and the measure of quadratures covariance. In 
both cases we evaluate the lower bound to the uncertainty in the phase-covariance-estimation given by photon noise, 
in function of the fundamental parameters: the interferometers central phases (z = 1,2), the quantum and 
classical beam intensities A and detection efficiency r]. The Results are reported in Sec. |V] For rather challenging 
conditions, namely almost ideal efficiency and perfect control of the stability of the interferometers central phases 
(regime (A)), TWB could deliver extraordinary advantage due to its photon number correlation at the quantum level 
(entanglement). This regime corresponds to the situation analyzed in [S]. Conversely, there exists a less demanding 
regime ((B)-regime described before), in which quantum strategies provide good enhancement in a more favorable 
experimental condition. In this case,both for TWB and SQB the expression of the minimal uncertainty presents the 
usual scaling with losses, (ex 1 — zy), and with the quantum light intensity, (oc 1/A), typical of single phase estimation 
using strong local oscillator and squeezed light. 

Finally we draw the conclusions in Sec. |VI| 

II. THE INTERFEROMETRIC SCHEME 

Let us consider a system as depicted in Fig[^ Two interferometers {i = 1,2) are injected at the ports denoted 
by the mode annihilation operators bi by a couple of identical coherent beams while the remaining ports 

identified by the mode operator Oi (unused in the classical scheme) are fed with a quantum state |4'(A))oi,02) where 
A is the mean number of photon in each mode. The readout ports are denoted by the mode operator Ci which will 
be function of the phases shifts (j)i among the arm of each interferometer, Ci = Therefore, a final combination 

of the outputs results in an observable C {ci,C 2 ,h.c.) = C ((/)i,(/) 2 )- A proper choice of the operator C leads to an 
estimation of the phase-noise correlation. Here, it is useful to recall the properties that the input-output operator 
relations of a linear interferometer (for example a Michelson-type) are equivalent to the ones of a beam splitter (BS) 
with transmission coefficient r = cos^(^/2). 

The losses in the system are taken into account by considering in both channels two identical detectors with quantum 
efficiency zy, formally re-defining the output operators with the substitution Ci —> y/fjci in the normal-ordering products. 
For example the photon number operator will be = r]c\ci. 

If only the classical field is injected, the photon counting statistics at the output ports is simply the one of coherent 





FIG. 1: Ideal scheme of the proposed experiment. The two modes of a twin beam are mixed with two identical coherent states 
at the input beam-splitters of two coupled interferometers (bolometer). The covariance of the phase noise at the outputs is 
studied in order to detect very faint (but correlated) noise. 


beams after reflection probability (1 — r^) and detection probability 77 , 

= va-nh, ( 1 ) 

{SNXt = {NXt = wii-n), 

{SN.SNXtr, = 0 - 

Now we consider two possible quantum states feeding the free input ports of 2i and two related readout strategies. 
Readout strategy 1: TWB ~ The TWB correlated state can be expressed in the Fock bases {|m)a.} as 


|m,TO)ai,a 2 ( 2 ) 

where 0 is the global phase, which we set in the following to 0 = 0 without loss of generality. The TWB presents 
perfect correlations in the photon number rrii = a\ai meaning that 01 , 02 ((TWB|(7Va^ — iVo 2 )'^|TWB))ai,a 2 = 0, 
VM > 0. It implies for example that variance of the photons numbers difference {S{mi — 777 , 2 )^) (with Sm = m — {m)) 
is identically null if losses are neglected. It also suggests to choose the measurement operator in the same form 
= {{Ni — N 2 X), since this should correspond to a reduction of the photon noise in the measurement, 
Anally improving the sensitivity. 

By using the equivalence between interferometers and BSs mentioned before, we can calculate the photon statistics 
of TWB transmitted to the output ports (in absence of classical coherent field) and detected with quantum efficiency 
77 . The mean photon number, the variance and the covariance are respectively: 

= gnX (3) 

= VTiX{l + r]TiX) 

{5N^5N2)™X = 0'Tir2A(I + A) 

Readout strategy 2: two squeezed states - The product of two single mode squeezed vacuum states writes: 

lOai 0 I0a2 = ^ai (05^2(0 10)ai ® |0)ii2 

where = exp[^ (“l)^ ~ is the squeezing operator. If we set ^ then A = sinh^ |^| represents 

the average number of photons of the squeezed vacuum, taken equal in both the modes. 

Defining the the quadrature of the field as Xi = and yi = , and supposing yi the squeezed and Xi the 

anti-squeezed one, it is known that in the single interferometer the injection of the squeezed field provides a fixed 
factor {5y1) = of resolution enhancement for arbitrary brightness of the coherent beam [T]. It is expected that 

the increased resolution in the estimation of the phase shifts 4>i ^md (/)2 separately reflects in a better estimation 
of their correlation if the correlation of the squeezed quadrature W of the output modes are considered, namely 
C = Xi- X2. 


\'^W)c 
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III. CORRELATIONS AT THE READOUT PORTS 

As a figure of merit for the correlations at the read-out ports, we study the noise reduction parameter NRF± = 
{d{Ni ± N 2 )'^)/{Ni + N 2 ) [TTMl4j . It represents the ratio between the variance of the photon number sum (difference) 
and the corresponding shot noise limit. The NRF- < 1 it is a well known condition of non-classicality for the 
correlations of a bipartite state and its value also determines the quantum enhancement achievable in certain sensing 
and imaging protocols [6]. For the same reason, NRF+ < 1 could be be interpreted as strong signature of anti¬ 
correlation of the photon number beyond classical limits. The BS-like transformation allows evaluating the fluctuation 
of the fields at the output ports in function of the input field. In particular one gets 

(TV.) = (iV)™ + (iV)“^ (4) 

(SNf) = (<5iV2)™ + (<5Ar2)“^ + 2(iV)™(iV)“^ (5) 


and the explicit expression, function of the parameters, can be directly obtained by substituting the quantities 
according to Eq.s Q and (|^. In particular Eq. (|^ shows that the covariance is composed by the TWB covariance 
{dNidN 2 }'^^^ and a second term containing the phase of the coherent field ijj, originated by the BS interaction of 
the two fields. Interestingly, the choice of ^ = 7r/2 maximize the covariance, while for '0 = 0 the covariance can even 
be negative (anti-correlation of photon numbers). _ 

The NRF± can be easily calculated from Eq.s ( |4|5|6[ ) by exploiting the identity {S{Ni ± N 2 Y) = {5Nl) + ((5A|) ± 
2{5Ni5N2)■ Using the same notation we have 


{6{Ni ± N2f)'^^^ 

NRF± - - 


2(iV)0f (1 + 2(iV)™ ± 2^/(^ivi^iv2)™ cos[20]) 
2(A)“'* + 2{N)TWB 


( 7 ) 


For simplicity, in Eq. Q and in the following we have assume ti = T2 = t. We note that in general the NRF- 
is minimized for tp = irj2 (corresponding to the optimization of the photon number correlation), while NRF^ is 
minimized when '0 = 0 (corresponding to the optimization of the anti-correlation). The NRF-{ip = 7r/2) and 
NRF+{'ijj = 0) are plotted in Figj^ In order to analyze the behaviour shown in the figures, and for the forthcoming 
discussion of the results concerning the phase-covariance estimation in Sec. |Vj it is useful to distinguish two regimes: 


Regime-(A): TWB-like correlations 


Referring to Eq.Q, when the intensity at the read-out port is dominated by the TWB, i.e. 
or explicitly k = /i(l — r) /rA ^ 1, the noise reduction factor reduces to the one of TWB, NRF± ~ {d{Ni ± 
iV2)2)™/2(Ar)™: 


NRF- ~ (1 — tit) + t]t{1 -I- 2A — 2a/A(1 + X))k ~ (1 — tjt). 


( 8 ) 


NRF+ ~ 1 -h r?T(2A -k 1) (9) 

recovering the expression of noise reduction factor for TWB in presence of losses [S]. Of course the condition <C 

appears quite challenging to achieve in the relevant case of practical interest in which the coherent mode is 
largely populated. Larger is closer to unity have to be the equivalent-transitivity r of the interferometers. In Fig. 
[^it corresponds to the region r ~ 1, where NRF_ drops to zero, clearly behaving as a singular point. We observe 
that in the same regime, as indicated by Eq. (|^, the NRF+ (marking the anti-correlation) rapidly grows accordingly. 
This regime has been studied in |8]. Even if in principle it allows exploiting the perfect TWB-like correlation, with a 
large classical power circulating into the interferometer, thus obtaining surprising quantum enhancement in the phase 
correlation estimation, the fragility of this regime has demanded for a more extended exploration of the parameter 
space. This is one of the motivation of this extended work. 
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(b) 
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FIG. 2: Noise Reduction Factor versus the transmission coefficient r of the TWB modes at the read-out ports of the interfer¬ 
ometers. Different colored-lines stand for different mean number of photons of the TWB (A): a) is the noise reduction factor 
{NRF-) exploiting quantum correlation at the output ports, for ip = n/2, while b) is the noise reduction factor {NRF^) 
exploiting quantum anti-correlation, for 7/1 = 0 . Here we consider the ideal lossless case, 77 = 1, and we set /i = 10®. 


Regime- (B): bright quantum correlation 


When the coherent power reflected to the measuring port is much higher than the transmitted power of TWB 

{Nr^f »(fv)™, 


i.e /r (1 — r) >> Ar, Eq.Q reduces to 

NRF± ~ 1 + 2(Ai)™ ± 2^{5N^5N2)™b cos[2ijj]. 

Introducing the explicit expressions of the various moments of the photon number distribution distributions we have 


NRF_{'P = 7r/2) = NRF+{'P = 0) ~ 1 -b 2?7r (a - ^^(l + A)) ~ 1 - ryy -b 


rjT 

4A’ 


( 10 ) 


where in the last equality we have also considered the limit of high intensity TWB, i.e. A ^ 1. It is worth to be 
notice that the NRF, for the proper choice of the phase of the classical fields, is always smaller than 1, whatever 
the intensity of TWB and losses. Thus, Ni and N 2 are always correlated (or anti-correlated) beyond the classical 
limit. It is possible to switch between quantum correlation and quantum anti-correlation just by acting on the phase 
Ip of the classical fields. Even more interesting, the correlation can be extremely bright, because the mean number 
of photon at the read-out ports is determined by the brightness of the classical beam = 77(1 — t)^, which can 

be increased almost arbitrarily in real experiments. It is clear from Eq. (101 that the highest correlation is obtained 
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when A ^ 1 and at the same time r ^ 1. For example in Fig. [^for plausible values of r ^ 0.90 and A = 10 the NRF 
is 0 . 1 , and the mean intensity of the output signal is (1 — T)/i = 10 ^ photons. 


IV. ESTIMATION OF PHASE CORRELATION (HOLOGRAPHIC NOISE) 

Since we the phase fluctuations due to the holographic noise are expected to be extremely small, we can expand 
C{ 4 >i, 4 > 2 ) around the chosen central values 4 ' 2 , 0 j namely: 

1 . 


Ci4’l^4’2) — Ci4‘l,OT4‘2,o) + ^i 0.C(0i_Oi ^2.o) 

6cl>^6<P2+0{d<l>^) 


( 11 ) 


where S(j)i = 4>i — 4>i o, and d^t\kC{(t)i o ,</>2 o) is the ih + fc)-th order derivative of C{(pi, 4 > 2 ) calculated at (j^i = 4>i O) 

hj = 1,2 

In order to reveal the HN, the bolometer exploits two different configurations: the one, “||”, where HN correlates 
the interferometers, the other, “_L ”, where the effect of HN vanishes. The statistical properties of the phase-shift (PS) 
fluctuations due to the HN may be described by the joint probability density functions /y (</>!, </> 2 ) and /_l(</>i, t!> 2 )- We 

make two reasonable hypotheses about fx{.4’ii4>2), x =||,-L. First, the marginals !Fx\4>i) = / <Pj)j hj = 1,2 

with i 7 ^ j, are exactly the same in the two configurations, i.e. one cannot distinguish between 

the two configurations just by addressing one interferometer. Second, only in configuration “_L” it is f±{4>i-, (f’ 2 ) = 
J-^\4’i)J^±\4>2), i.e., there is no correlation between the PSs due to the HN. Now, the expectation of any operator 
0 ((/)i,(/) 2 ) should be averaged over fx, namely, (/) 2 )) £x 0 (</>i,()) 2 ) = /(0((/)i, (/) 2 )) fx{(t>i,(t> 2 ) d^i d(j) 2 . In 


turn, by averaging the expectation of Eq. (|I1|), we have: 


C{4>1,4>2) — (^(01,0, ^2,0)) + 2^* (^0i,0i^('('l,O, <('2.0)) ^x 

+ {dl,^,C{<Pi,o,h,o)) £x [Scj,,6cl>2]+0{6cl>^) 


( 12 ) 


where we used £ 3 , [Scj)^] = 0. Then, according to the assumption on fxi(l>i,4’2) we have = £_l \p4>V\ 

f_L [5(j)i5(j)2] = [54'i\£i. [()</> 2 ] = 0, and from Eq. (12) follows that the phase-covariance may be written as: 


£\\ [54>i5(j>2] 


C{ct>iA2) 




C( 01 , 02 ) 


(^^1.02^('('1,O, 4>2fi)) 


(13) 


that is proportional to the difference between the mean values of the operator C( 0 i,(^ 2 ) as measured in the two 
configurations ‘jj and “T”. 

Indeed, one has to reduce as much as possible the uncertainty associated with its measurement: 


U{5(j)i54>2) 


Vary 

C(<^ 1 ,<^ 2 )' 

-1- Var_L 

'C{^i,<l>2) 


i 


where Var^ 


C'(</)l,</. 2 ) 


= f 


(</'!, <(' 2 ) 


_ f 


(^0i,02^(‘(’1>O’*(’2.o)) 


C((/<l,</- 2 ) 


^ 2 


{54>i,64>2 < 1 ) 


(14) 


Under the same hypotheses used for deriving Eq. (13) we can calculate the variance of < 7 (^ 1 , ^ 2 ) as 

'C{ct>^,ct>2) 


Var,, 


= Var 


</' 2 ,o) + Sfe Akk £x + ^12 £x [^<Pl5(t>2\ + 0{S(j)^) 


(15) 


where: 


Akk = (C'(^i,o, <)>2,o)i90fc,0fcU(0i,o, (>^ 2 , 0 )) 

+ {[d^k^i4'i,o,4‘2,o)]'^) - {C{4'l,0,4‘2,o)){dl^^^^C{(j)i^o,4>2,o)) 

Ai2 = 2 (C'((/)qo, ^(>2,0)3^1,02^(^i.o, ^2,0)) 

-I- 2 (i 90 ^C'(())i,o, <)'2,0)1902C'(</^i,o,</^’2.o)) — {C{4>i,o,(l>2,o)){dlj^^^^C{(j)i^o,(j)2fi)) 


(16) 
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Analyzing expression (151, we note the presence of a zeroth-order contribution that does not depend on the PSs in¬ 
trinsic fluctuations, and represents the quantum photon noise of the measurement described by the operator C{(pi,(j) 2 ) 
evaluated on the optical quantum states sent into the bolometer. The statistical characteristics of the phase noise enter 
as second-order contributions in Eq. (15) from each interferometer plus a contribution coming from phase correlation 
between them. 

This work addresses specifically the problem of reducing the photon noise below the shot noise in the measurement 
of the HN, therefore in the following, we will assume the zero-order contribution being the dominant one. Of course, 
this means to look for the HN in a region of the noise spectrum that is shot-noise limited. Since the HN is expected 
up to frequencies of tens MHz, it follows that all the sources of mechanical vibration noise are suppressed. Therefore, 
the zero-order uncertainty that we will study here is 




2 Var 




(17) 


A. TWB 


As we argue in Sec. [Hj when TWB state is injected it should be promising to define the observable operator in 
the form C {(j)i,4>2) = (IVi — A^ 2 )^ , (Af > 0) because of the perfect photon number correlation of TWB. Indeed, in 
the Sec. |HI| we show that at least up to the second order M = 2, the strong non-classical correlations are preserved 
at the output ports of the interferometers (for ip = 7r/2), justifying the conjecture that an advantage in terms of 
noise reduction would be obtained if we can estimate the phases covariance starting from the measurement of an 
observable of that form. We notice immediately that for M = 1, corresponding to the photon numbers subtraction, 
the proportional coefficient in Eq. (131, containing the double derivative with respect to both the phases will be null. 


Thus, we have to move to the second order measurement i.e. C{4>i,(j)2) = (A^i(</i’i) — A^ 2 (<(> 2 ))^ = Nf + N 2 — 2 N 1 N 2 . 
Hereinafter we also consider the same central phase of the two interferometers ^ 1.0 = <(' 2,0 = 4>o- 
According to Eq. (IT^ we get: 


£\\ [6(j)id(l)2] 


[iViiV2] - £^ [iViiVa] 


(18) 


where we have used again the symmetry of the statistical properties of the two interferometers, in particular 
'S’ll(-L) [-^ 1 ] = ^11 (_L) [-^ 2 ]- The covariance of the phase noise is proportional to the difference between the photon 
number correlation when the phase noise is correlated (||) and when it is not (T), as one could expect. 

The uncertainty of the measurement, due to photon noise can be obtained by Eq. (17) where \&r\Ci(bi n, (1)2 o)l = 

m{M-N2{<i)o)f)-m{M-N2{(i)o))Y- 


B. Independent Squeezed States 


It is rather intuitive that the most simple form of the measurement operator C{(j)i, ^ 2 ), that combines the squeezed 
quadratures measured at the read-out port, and has non-null mixed derivative with respect to the phases, ^0, 

would be the product hi -12 (where Yi are the squeezed quadratures). However, to avoid the presence of a dc-component 
in the measurement it turns out more useful to consider the fluctuation of the quadratures around their central value, 
therefore defining C = (yi((()i)—£[yi])-(F 2 (<(' 2 )—'£’[^ 2 ]), where we have taken into account that £\\ \Yi\ = = £\Yi]. 

The covariance of the phases is estimated according to Eq. (131 as: 


£\\ [5(j)i5(j)2] 


£\\ [Y1Y2] - £± [^1^2] 

{dl,^YiiMY2{(l)o))’ 


(19) 


Since the fluctuations of the quadratures due to quantum noise are independent in the two interferom¬ 
eters, the zero order uncerainty on the measured observable (ses Eq. (17)) remains Var[C((/)i_ 0 ; </' 2 ,o)] = 

m{M-£[Y,]f){iY2i(l)o)-£[Y2]f). 

















V. RESULTS 


The calculation of the variance of the measurement operator C(4io), in particular for TWB case, involves many 
fourth-order terms of the photon number operator, i.e. eighth-order product of field operator Ci and c|, and the 
calculation and the complete expression of this variance are too cumbersome to be reported here. Thus, we will 
present numerical results for the most significative regions inside the parameter space and we give some general 
expression in particular relevant limits. First of all we need to define the classical benchmark to compare performance 
using quantum light. The uncertainty achievable in the estimation of the phase covariance, if only the coherent beams 
are used, is = y/2/{r]^cos'^ [4 >q/2]). We notice that the scaling is with the detected number of photons, i.e. the 
square of the shot noise limit typical of the single phase estimation. This is because we are measuring a second 
order quantity, namely the covariance of the phases. As usual, it is clear that without any particularly low energy 
constraint, in order to reach high sensitivity in a phase-correlation measurement is necessary to push the intensity of 
the classical field. Therefore, even the quantum strategy should face and should improve the sensitivity when high 
power is circulating into the interferometers. In general we will consider the limit ^ ^ 1. 

Concerning the use of the two independent squeezed states, we can summarize the results in the following couple 
of equation; 


^SQ /^CL 


2 4A > 

1 - ry(l -I- cos[(/)o])a/A(1 - a/A) (/r>l,A<Cl) 


( 20 ) 


Of course, we expect that the advantages of using squeezing, and in general quantum light, is effective when it 
experience a low loss level. Thus the most interesting regime is when the two interferometers transmit almost all 
the quantum light to the read-out port, meaning that the ce ntra l phase must be close to 0, according to the BS-like 
behaviour r = cos^[())o/2] — 1. Applying this limit to Eq.s (20 


we have 


( 0 ) 

CL 


I — T] + ?7/4A for a 1 and 


^SqV^cl ~ 1 ~ — y/X) for A <C 1. We can appreciate visually what does it mean by looking to the Fig. 

A flat region (in logarithmic scale) appears in the uncertainty reduction in function of the central phase ^o- Since 
A = 10 in the hgure, the value of the uncertainty reduction, given by the previous expression, is well represented by 
1 — ?? -b 77 / 4 A. In the opposite limit of A ^ 1 the advantage of squeezing is lost, according to the second of Eq.s (20|. 
For example for A = 0.1, represented in hgure Fig. the improvement is reduced to 0.6. 
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FIG. 3: Uncertainty provided by the twin beam state (TWB) and by the product of squeezed states (SQ) normalized to the 
classical limit in function of the central phase (f>o in which the interferometers are operated. The plot is obtained for A = 10, 
H — 3 X 10^^, tjj = 7r/2, and different colors corresponds to different values of the detection efficiency 77 . 

Concerning TWB, one can clearly discern two different regions both in Fig. [^and Fig. one for really small values 
of the central phase, namely (j)Q < 10“® and an other one in the range 10“® < (po < 10 . They correspond, for the 

a specihc choice of the parameter indicated in hgure, to the two relevant regimes that have been individuated in Sec. 
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FIG. 4: Uncertainty provided by the twin beam state (TWB) and by the product of squeezed states (SQ) normalized to the 
classical limit in function of the of the central phase (po in which the interferometers are operated. The plot is obtained for 
T] = 0.95, fi = 3 X 10^^, Ip = 77 / 2 , and different colors corresponds to different values of the mean number of photons per mode 
A of the quantum light. 





FIG. 5: Uncertainty provided by the twin beam state (TWB) and by the product of squeezed states (SQ) normalized to the 
classical limit in function of the quantum efficiency t]. The plot is obtained for (p = 10“®, /r = 3 x 10^^, ip = 7r/2, and different 
colors correspond to different values of the mean number of photons per mode A of the quantum light. 


(A)-TWB like correlations- when <C . The condition is guaranteed if the central phases are close 

enough to zero, (pi^o = (p 2 fi — Oj meaning that the transmissivity of the equivalent-BS approaches the unity 
Ti = Cos^[(pi^Q/2] ~ l(t = 1,2). This is the regime studied and reported in [5]. For intense coherent beam and 

intense TWB source, i.e. ^ A 3> 1, one gets ^twb/^ctI ~ 2-\/5(l — rf), while in the case of faint TWB, 
A <C 1 and ^ » 1, the result is ~ 1 / 2(1 — ? 7 )/? 7 . In both cases TWB allows reaching an amazing 

uncertainty reduction that approaches zero for unitary detection efficiency. This behaviour is clearly shown in 
Fig. [5] in which the choice of /x and (pQ ensure to be in the TWB-like regime at least for the range of values 
of A represented there. All the TWB curves drop to zero, and for some value of the efficiency, depending on 
the intensity A, they fall below the corresponding SQ curves. The limits for A 1 and for A » 1 are also 
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reported in dotted and dashed black lines respectively. However, we observe that for quantum light intensity 
A > 1 reachable in experiments nowadays (for example A = 3 in the picture) squeezing performs far better than 
TWB except for extremely demanding overall detection efficiency. 

(B)-Bright quantum correlations- when ^ . This regime corresponds to the flat region shown in 

Fig.sH] and 1^ for intermediate values of the central phase (f)Q. Aside a constant factor, the uncertainty reduction 
ior 1 behaves as for the two independent squeezing case, specifically 

easily appreciated when comparing the corresponding curves for TWB and Squeezing in the figures (taking in 
to account the logarithmic scale). 

For the sake of completeness, we did the same analysis considering to exploit the anti-correlations, dehning the 
observable as C(0i,(/)2) = (A^i -I- A^ 2 )^ = A^i + A(f + 2 N 1 N 2 (for ijj = 0) instead of the correlation when TWB are 
injected, obtaining analogous results in the regime of /r (1 — r) >> At. 

VI. DISCUSSION AND CONCLUSIONS 

In Sec. |V]we observed interesting features, leading to promising experimental conditions. Referring to Fig. [^and 
1^ there is an extended range of value of the central working phase (/>o of the interferometers in which the uncertainty 
reduction achievable by adopting quantum light is stable, at the value ~ 1 — ^ + ^/4A both with SQB and 

TWB (a part a factor in the last case). This kind of scaling is a well known results of phase estimation in a single 
interferometer combining coherent strong field and single mode squeezed light (in fact 4A « in the limit 4A ^ 1) 
[T]. Therefore, it turns out that a measurement of the phase correlation retains the same advantage of the single 
phase estimation. As an example, for rj = 0.9 and A = 3, compatible with the actual technology, we have 5.7 times 
of uncertainty reduction in the single measurement. Since in any experiments Af measurements are performed and 
the final uncertainty is U/VN, one would easily obtain the same sensitivity with a number of runs 30 times smaller, 
hence reducing the total measurement time of the same amount. 

On the other side, only for TWB, there exists a special setting of the central phase of the interferometer, when 
the classical fields component is made negligible at the read-out ports with respect to the TWB component, in 
which in principle the uncertainty reach the zero point, whatever the intensity of TWB. In particular for faint 
TWB (A <C 1) the uncertainty scales as Z-Z^wb/^cl ~ \/2(l — ifj/r] while for intense TWB (A 1) one gets 

^twb/^CL ~ (1 — ??)■ Even if at first glance this looks rather exciting, Fig. shows that in terms of absolute 

sensitivity squeezing performance can be overtaken only for rather high detection efficiency. For example for A = 3 
we expect Z^4wb ^SQ ^ — 0-99- However, for limited quantum resources, namely A < 1, TWB performs better 
than squeezing already for smaller and more realistic efficiency values. 

In conclusion, we have analyzed in detail a system of two interferometers aimed at the detection of extremely faint 
phase-fluctuation. The idea behind is that a correlated phase-signal like the one introduced by the “holographic 
noise” could emerge by correlating the output ports of the interferometers, even when in the single interferometer it 
confounds with the background. We demonstrated that injecting quantum light in the free ports of the interferometers 
can reduce the photon noise of the system beyond the shot-noise, enhancing the resolution in the phase-correlation 
estimation. Our results basically confirms the benefit of using squeezed beams together with strong coherent beams in 
interferometry, even in this correlated case. However, mainly we concentrated on the possible use of TWB, discovering 
interesting and probably unexplored areas of application of bipartite entanglement and in particular the possibility 
of reaching in principle surprising uncertainty reduction. 
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